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ON  THE  RELATION  BETWEEN  THE  VELOCITY  COEFFICIENT 
AND  BOUNDARY  VALUE  FOR  SOLUTIONS  OF  THE 
ONE-DIMENSIONAL  WAVE  EQUATION  * 

R.  MICHAEL  LEWIS  AND  WILLIAM  W.  SYMES  1 

Abstract.  The  one-dimensional  acoustic  wave  equation  is  a  simple  model  of  wave  propagation 
in  layered  media.  As  such,  it  is  used  in  theoretical  seismology  to  study  the  relation  between  sound 
velocity  and  the  surface  response  or  seismogram  (boundary  value)  as  a  simple  instance  of  the  reflection 
seismology  problem. 

Questions  about  the  nature  of  the  dependence  of  boundary  values  on  the  velocity  coefficient  in  the 
wave  equation  arise  naturally  in  this  context,  particularly  in  connection  with  perturbational  techniques. 
We  study  the  dependence  on  the  velocity  coefficient  of  boundary  values  for  solutions  of  the  one¬ 
dimensional  wave  equation.  We  present  estimates  for  both  the  Lipschitz  continuity  and  the  linearization 
error  for  the  map  between  velocity  coefficient  and  boundary  value.  In  particular,  we  show  that  this 
relation  is  Lipschitz  continuous  for  velocities  in  H2  and  differentiable  for  velocities  in  H3 .  We  also 
discuss  the  anomalous  smoothness  of  this  map  in  “oscillatory”  directions,  which  helps  explain  a  key 
idea  in  reflection  seismology. 

Key  Words,  wave  propagation,  one-dimensional  acoustics,  inverse  problems,  stability 
AMS(MOS)  subject  classification.  Primary  35L05,  secondary  35R30,  73D99,  86A15 

1.  Introduction.  Various  wave  propagation  problems  in  a  plane-stratified  half¬ 
space  {z  >  0}  can  be  modeled  by  the  one-dimensional  acoustic  wave  equation: 

(1)  -  alu  =  o. 

In  this  model  the  mechanical  properties  of  the  medium  are  described  by  the  sound 
velocity  c(z). 

We  impose  the  boundary  and  initial  conditions 

(2)  dzu(0,t)  =  f(t) 

(3)  u(z,t)  =  0  t  <  0, 

where  the  Neumann  datum  /  vanishes  for  |f|  large  so  that  the  disturbance  is  transient. 
The  velocity  coefficient  c  is  a  strictly  positive  function  and  for  the  moment  is  assumed 
to  be  C°°. 

If  the  datum  /  is  C°°  then  we  can  sensibly  define  a  map  Ff  :  C°°  — ►  C°°  by 

(4)  Ff  :  c  i->  dtu(0,t). 

If  u  represents  displacement,  then  the  map  Ff  relates  the  velocity  coefficient  and  the 
surface  velocity.  In  the  context  of  the  reflection  seismology  method  in  oil  exploration, 

*  Research  sponsored  by  AFOSR-89-0363.  A  version  of  this  work  appeared  in  R.  M.  Lewis’  doctoral 
thesis,  written  under  the  supervision  of  W.  W.  Symes  in  the  Department  of  Mathematical  Sciences, 
Rice  University. 

*  Department  of  Mathematical  Sciences,  Rice  University,  Houston,  Texas,  77251-1892. 
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Ff  represents  the  relation  between  the  surface  observations  (“seismogram”)  and  the 
unknown  subterranean  material  properties. 

The  conditions  (1),  (2),  and  (3)  form  a  well-posed  initial-boundary- value  problem  in 
the  half-space  {z  >  0}  for  the  distribution  datum  /  =  —  6  since  the  velocity  coefficient  is 
sufficiently  smooth.  However,  in  this  case  we  qualify  the  definition  of  the  map  between 
velocity  and  boundary  value  given  by  (4)  to  account  for  the  absence  of  smoothness  of 
/  =  —8.  The  progressing  wave  expansion  for  u  (see  Section  2.3)  shows  that  dtu  has  a 
well-defined  trace 


dtu( 0,  t )  =  constant  •  8{t )  +  v(t), 
where  v  €  C°°[0,  T]  for  any  T  >  0.  We  define  the  map  F  to  be 
(5)  F  :  c  t-+  u(t). 

Since  for  smoother  choices  of  /  we  have 


Ff(c)  =  f  *  (constant  •  8  +  F(c)) , 


properties  for  Ff  follow  from  those  of  F  via  well-understood  properties  of  convolution. 
Accordingly,  we  restrict  attention  to  the  case  /  =  —8,  the  so-called  impulse  response 
problem. 

To  this  point  we  have  assumed  that  the  velocity  coefficient  c  is  C°°.  Our  goal  in 
this  paper  is  to  establish  weaker  regularity  conditions  on  c  under  which  F  extends  to  a 
continuous  map  and  a  differentiable  map.  Our  main  results  are  that 

(a)  F  is  Lipschitz  continuous  as  a  map  H2  — >  L 2,  and 

(b)  F  is  differentiable  as  a  map  H3  — >  L 2  and  the  derivative  is  Lipschitz  continuous, 
as  we  shall  now  describe. 

Let  u  be  the  solution  of 


(6) 


dttu  -  d2zu 

=  o, 

z  >  0 

<92u(0,  t) 

=  -m 

u 

=  o, 

t  <  0. 

A  consequence  of  the  estimates  we  will  present  in  Section  2.5  is  that  the  map  F  relating 
sound  velocity  and  boundary  value  via  (5)  extends  to  a  map  H1  —>  L2. 

Our  first  major  result  is  that  F  extends  to  a  Lipschitz  continuous  map  from  H 2  to 

L2. 

THEOREM  1.1.  Given  a  velocity  c  and  T  >  0,  there  exist  D  and  r  >  0  such  that  if 
8c  e  H2[0,  D ],  ||  8c  ||#2[0iD]  <  r,  then 

||  F{c  +  Sc)  —  F(c)  ||i2[0i-r]  <  A  ||  8c  ||//i [0,z?]  • 

The  constant  K  depends  on  the  depth  D,  ||  c\\H2y0Dy 

cmin  =  inf  c(z), 

0 <z<D  v  ' 
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and  r. 

Our  second  theorem  concerns  the  differentiability  of  F.  Formally,  the  first-order 
effect  of  changing  the  velocity  from  c  to  c+  8c  is  to  perturb  the  solution  u  by  Su,  where 
8u  is  the  solution  of  the  linearized  problem 


(7) 


~^d2tt8u  -  d2zz8u 
c2(z ) 

dz8u(0,  t) 

8u 


2  8c(z) 
(?(z) 


0 


0, 


z  >  0 


t  <  0. 


As  with  the  trace  dtu( 0,  t )  there  is  a  highly  singular  part  and  a  smooth  part  to  dtSu( 0,  t). 
The  progressing  wave  expansion  for  Su  shows  that 

dtu( 0,  t )  =  constant  •  8\t)  +  constant  ■  8(t)  +  v(t), 

We  will  consider  a  formal  linearization  DF  of  F  defined  by 


DF[c]  :  8c  ►  u(0,f). 

In  order  to  determine  whether  DF  really  is  the  derivative  of  F,  we  will  estimate  the 
Taylor’s  series  remainder 


F(c  +  Sc)  -  F(c)  -  DF[c\(8c) 


for  all  velocity  perturbations  8c  which  are  “sufficiently  small.” 

Our  second  result  —  the  more  significant  for  applications  —  is  that  F  is  differen¬ 
tiable  as  a  map  H3  — >  L 2 ,  and  the  derivative  is  Lipschitz  continuous. 

Theorem  1.2.  Given  a  velocity  c  and  T  >  0,  there  exists  D  such  that  F  is 
differentiable  at  c  as  a  map  F  :  H3[0,D]  — ►  L2[0,T],  and  the  derivative  is  Lipschitz 
continuous. 

More  precisely,  given  c  and  T  >  0,  there  exist  D  and  r  >  0  such  that: 

(i)  For  all  8c  €  H2[0,D],  ||  8c  ||#2[0  £>]  <  r,  the  linearization  error  satisfies 

||  F(c  +  Sc)  —  F(c)  —  DF[c\(8c)  ||L2[0,t]  —  K  I!  IIl2[o,p]  II  ll/f2[o,£>]  F  B\\8c  ||i2[0jD]  • 
The  constant  K  depends  on  the  depth  D,  ||  c  ||^2[0  Dy 


C min  —  inf  c(z), 

o/ ^ / n  '  ' 


0  <z<D 


and  r,  while  the  constant  B  depends  on  the  depth  D,  ||  c  || ^3[0  and  r. 

(ii)  The  derivative  depends  continuously  on  c  in  the  sense  that  if\\  c  —  c  ||#2[0£>]  <  r, 

then 


II  DF[c]  -  DF[c\  ||  <  B\\c-c  \\H2[0M  <  r, 

where  the  constant  B  depends  on  the  depth  D,  ||  c||//3jo£>p  cmin>  and  r,  and  the  operator 
norm  is  taken  over  operators  from  H2[0,D ]  to  T2[0,T]. 
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The  appearance  of  ||  c||#3[0D]  in  the  regularity  estimates  suggests  that  F  is  differ¬ 
entiable  with  Lipschitz  continuous  derivative  only  as  a  map  H3[0,D ]  — >  L2[0,T],  The 
smoothness  restriction  c  €  H3  is  a  consequence  of  the  singular  nature  of  the  boundary 
datum  /  =  —8.  For  smoother  /  we  can  relax  the  smoothness  requirement  on  the  ve¬ 
locity  at  the  expense  of  the  appearance  of  derivatives  of  /  in  our  estimates.  For  such  a 
resale,  see  [1]. 

The  particular  form  of  the  linearization  error  estimate  has  further  significance.  On 
the  one  hand,  the  error  estimate  indicates  that  linearization  is  valid  only  for  velocity 
perturbations  with  small  if2-norm.  On  the  other  hand,  it  turns  out  that  we  can  be 
much  more  optimistic  about  the  accuracy  of  the  linearization  of  F  in  the  special  case 
where  8c  is  highly  oscillatory.  The  key  point  is  that  the  map  F  is  very  smooth  in 
these  directions ,  as  we  shall  now  discuss.  This  observation  is  significant  for  the  use  of 
the  problem  (6)  as  a  model  in  reflection  seismology  and  helps  explain  an  item  of  the 
conventional  wisdom  in  that  subject. 

To  simplify  notation,  let  ||  •  ||fc  denote  the  norm  on  Hk[0,D],  with  H°[0,D]  = 
T2[0,D].  Suppose  that  the  support  of  8c  is  contained  in  (0,D).  We  will  say  that 
8c  is  oscillatory  if 

f  8c  <||  6c||0. 

Jo  0 

For  convenience  we  will  write 

/  6cl  =£  Iloilo. 

Jo  Ho 

where  e  should  be  thought  of  as  a  small  number.  For  our  discussion,  we  will  also  assume 
that  the  H3- norm  of  8c  satisfies  ||  8c  ||3  <  R. 

If  8c  is  oscillatory,  then 

||<5c||q  =  J  dz  8c  (z)  (J  8cj  <  |  J  8c 
and  consequently 


II  Mo  <eW\i- 

Thus  if  8c  is  oscillatory  we  have 

II Mo || M2  <  IIM|0(||M|1  +  Iloilo) 

<  II  llo  (ll  Hi  +  II  Sd  llo  II  f>c'"  Ho) 

<  II M0 II M?  ||  <5c  ||| 

<  ||  8c  ||i  . 

In  this  case,  the  estimate  of  the  linearization  error  in  Theorem  1.2  reduces  to 

II  F(c  +  «c)  -  F(c)  -  DF{c](6c)  ||1![0iT]  <  CeRl  ||  Sc ||f . 
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The  point  is  that  the  error  estimate  is  superlinear  in  the  Ft1- norm  of  5c  as  long  as  8c 
is  oscillatory.  Moreover,  the  oscillation  of  6c  enables  us  to  control  the  size  of  the  term 

CeRK 

As  an  example  of  a  family  of  oscillatory  perturbations,  let  x  €  Qj°(0,  D),  0  <  x  <  1, 
and  x  =  1  on  some  interval  [20,  C  (0,  D ).  For  e  >  0  and  a  >  1  define 

6ce(z)  =  £ax{z )  sin  -• 

£ 

Then 

=  0(e)  ||  6ce  ||0 
0 

while 

II  lit  =  0(ea~k). 

In  particular,  if  a  >  1  then  6c£  G  H1,  and  if  a  >  1  then  ||  6ce  ||x  — >  0  as  e  — *•  0. 

If  a  =  3,  6cc  stays  in  a  bounded  set  in  H 3  as  e  — ►  0.  On  the  other  hand,  for 
1  <  o.  <  3,  ||  8cs  ||3  — >  00  as  e  — »  0.  Nonetheless,  since  6ce  G  H1,  F(c  +  ScE)  is  still 
defined  and  we  have  the  estimate 

||  F(c  +  Sc)  -  F(c)  -  UF|c|(&)  ||1![0  T]  <  Ce ^  ||  Sc  ||f  . 

Since  a  >  1,  the  quantity  Ce~5~  —■»  0  with  e.  Thus  the  linearization  error  is  small,  even 
though  the  if3-norm  (and  the  H2- norm,  for  1  <  a  <  2)  of  6ce  is  unbounded  as  e  — ►  0! 
This  example  illustrates  the  anomalous  smoothness  of  F  in  the  directions  of  oscillatory 
perturbations. 

This  observation  helps  to  make  precise  an  idea  central  to  the  mathematical  methods 
of  reflection  seismology.  One  of  the  tenets  of  reflection  seismology  methods  in  petroleum 
exploration  is  that  linearization  of  the  map  F  (or  its  equivalent  in  more  complex  models) 
is  valid  as  a  means  of  studying  the  effects  of  oscillatory  perturbations  of  a  smooth 
reference  medium.  This  technique  is  the  basis  of  many  of  the  methods  used  to  analyze 
reflection  data  in  oil  exploration.  The  preceding  analysis  of  the  linearization  error  in  the 
case  of  oscillatory  perturbations  supports  this  belief,  and  gives  a  rigorous  justification  for 
the  “velocity/reflectivity”  dichotomy  that  makes  seismic  data  interpretation  possible. 

Another  important  property  of  the  map  F  is  the  effect  of  slowly-varying  perturba¬ 
tions,  the  obverse  of  the  case  we  just  considered.  The  effect  of  slowly-varying  pertur- 
baticnc  13  related  to  the  smoothness  restriction  c  G  H3  required  for  differentiability.  To 
understand  how  F  reacts  to  slowly- varying  perturbations  we  will  consider  the  closely 
related  impedance  problem.  The  issue  of  regularity  for  the  impedance  problem  was  stud¬ 
ied  by  Symes  in  [4].  There  he  considered  the  relation  between  coefficient  and  boundary 
value  for  Webster’s  horn  equation: 

(vix)9ft  -  dx(r}(x)dx))  v  =  0,  x  >  0 
dxv(0,t)  =  —6(t) 

=  0, 


v 


t  <  0. 
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The  quantity  7](x )  is  the  impedance.  If  we  set  cr(x)  =  c^Iog  r}{x),  then  we  obtain  the 
more  convenient  formulation 

(dtt  ~  dL  ~  <r(x)d*)  v  =  0,  x  >  0 

(8)  dxv(0,t)  =  —S(t) 

v  =  0,  t  <  0. 

In  this  case  the  coefficient  a  and  and  the  boundary  value  dtv(0 ,t)  are  related  by 
the  map  r: 


r  :  <7  i — ►  5tu(0,f). 

to  T2[0,T];  i.e.,  the  relation 
impedance  problem  via  the 


The  impedance  is  given  by 


Symes  showed  that  r  is  a  C1-diffeomorphism  from  L2[0,T] 
is  differentiable  for  impedances  in  H 1 . 

The  velocity  problem  (6)  can  be  transformed  into  the 
change  of  variable 


V(x)  =  c  x(z(;r)). 

This  change  of  variable  is  called  the  travel-time  transformation.  The  quantity  x  in  (9) 
is  called  the  travel-time ,  and  is  the  time  it  takes  a  signal  traveling  downward  from  the 
surface  to  arrive  at  depth  z. 

A  salient  difference  between  the  velocity  problem  and  the  impedance  problem  is 
the  effect  of  the  travel-time  transformation  in  “straightening  out”  the  characteristic 
curves  associated  with  the  differential  operators.  In  the  impedance  problem  (8)  the 
characteristics  are  independent  of  a.  In  the  velocity  problem  (6)  different  velocities 
give  rise  to  different  characteristics,  a  feature  that  enormously  complicates  the  proof  of 
the  regularity  theorem  and  forces  us  to  require  additional  smoothness  on  the  velocity  c  in 
order  that  F  be  differentiable.  This  simplifying  effect  of  the  travel-time  transformation 
(9)  reflects  a  fundamental  feature  of  the  map  F. 

As  previously  noted,  Symes  established  differentiability  on  H 1  for  the  impedance 
problem.  By  contrast,  for  the  velocity  problem  a  linearization  estimate  in  terms  of  the 
fifTnorm  of  6c  was  derived  only  for  oscillatory  perturbations  of  a  smooth  ( H 3)  velocity. 
That  no  better  result  is  possible  is  inherent  to  the  velocity  problem.  While  the  map  F 
is  surprisingly  well-behaved  for  oscillatory  perturbations  in  the  velocity,  it  is  extremely 
sensitive  to  slowly-varying  perturbations  in  the  velocity. 

The  extreme  sensitivity  of  F  to  changes  in  the  slowly- varying  component  of  the 
velocity  is  manifest  in  a  phenomenon  known  as  the  travel-time  shift.  The  travel-time 
shift  refers  to  the  possibly  catastrophic  effect  of  a  small  change  in  the  sound  velocity. 
The  arrival  times  of  reflected  signals  at  the  surface  are  determined  by  the  sound  velocity 
in  the  medium.  On  the  one  hand,  a  small,  slowly-varying  perturbation  in  the  sound 
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velocity  will  produce  only  small  changes  in  arrival  times.  However,  if  the  reflected  signals 
are  highly  oscillatory  this  can  lead  to  a  dramatic  change  in  the  surface  measurements, 
since  the  highly  oscillatory  features  of  the  signal  will  have  undergone  a  shift  relative 
to  their  former  locations.  The  additional  smoothness  requirements  of  our  regularity 
theorem  serve  to  ameliorate  this  unfortunate  behavior.  For  a  fuller  discussion  of  the 
travel-time  instability,  see  [3]  and  the  references  cited  there. 

The  remainder  of  this  paper  is  devoted  to  the  proof  of  Theorems  1.1  and  1.2.  The 
proof  is  rather  lengthy.  It  is  perhaps  possible  that  a  shorter  proof  might  be  possible 
via  a  study  of  the  travel-time  map.  However,  such  an  argument  would  be  inextricably 
bound  to  the  result  for  the  impedance  problem,  which  has  no  analog  for  problems  with 
multiple  wave  speeds  and  most  especially  for  the  multi-dimensional  wave  equation. 
This  consideration  motivates  the  proof  given  here  as  a  pattern  of  attack  on  the  multi¬ 
dimensional  problem  using  various  energy  estimates. 

The  proof  extends  the  methods  employed  by  Symes  in  his  study  of  Webster’s  horn 
equation  in  [4],  namely,  sideways  energy  estimates  and  the  progressing  wave  expansion. 
Sections  2.2  and  2.3  introduce  these  ideas.  The  application  of  these  techniques  to  the 
velocity  problem  (6)  is  complicated  by  the  more  complex  characteristic  geometry  of 
the  velocity  problem.  The  heart  of  the  proof  lies  in  Section  2.5,  where  we  use  the 
sideways  energy  identity  to  derive  trace  estimates  for  solutions  of  the  wave  equation. 
These  estimates  are  sharp  propagation  of  regularity  results.  In  Sections  2.6  and  2.7  we 
use  these  trace  estimates  to  compare  solutions  of  the  wave  equation  associated  with 
different  velocities.  The  remaining  sections  complete  the  proof  of  the  theorems. 

2.  Preliminaries  to  the  proofs  of  Theorems  1.1  and  1.2. 

2.x.  Time-depth  relationships.  In  this  section  we  will  we  will  define  the  depth 
D  and  distance  r  which  appear  in  Theorems  1.1  and  1.2  and  discuss  the  relationship 
between  the  time  and  depth  domains.  We  will  assume  that  F(c),  F(c+8c),  and  DF[c\(6c ) 
are  given  for  time  t  in  an  interval  [0,  T],  The  definition  of  a  related  interval  in  depth  z 
is  somewhat  complicated  —  here  we  encounter  the  first  of  many  difficulties  caused  by 
the  difference  in  the  characteristics  of  the  wave  operators  associated  with  the  different 
velocities. 

Given  c,  let  C\  —  c  and  c2  =  c  +  8c,  and  for  i  =  1,2,  define  the  two  travel-times 


Mz)  =  [  dC 

Jo 


1 

c.(0' 


Since  the  fa  are  increasing  functions,  they  have  well-defined  inverses. 
Given  a  time  T,  define  the  depth  d  by 


This  defines  a  ‘'domain  of  dependence”  for  the  surface  measurements  on  the  interval 
[0,  T],  If  signals  travel  with  velocity  c\,  then  d  is  the  deepest  point  a  signal  can  reach 
and  still  return  to  the  surface  before  time  T. 
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Now  choose  b  >  0  and  set 

c  =  ^(f  +  5)- 

With  b  and  D  thus  defined,  suppose  that  ||  c2  —  ca  ||fflj0Dj  is  small  enough  that 

||  <t>2  —  IIl°°[0,D]  <  2' 

Then  with  these  choices  we  have 

t  +  MD)  >  T  —  <f>\(D) 

r  +  UD)  >  T-MD) 


for  all  r  >  —b/ 2. 

The  depth  D  is  slightly  deeper  than  the  depth  reached  at  time  Tj 2  by  signals 
travelling  at  speed  c\.  The  slightly  enlarged  depth  domain  is  required  in  order  to 
compare  signals  travelling  at  a  slightly  slower  velocity  c2. 

Now  we  can  say  how  small  6c  must  be  in  order  for  the  estimates  in  Theorems  1.1 
and  1  0  to  hold.  Let  cmtn  and  cmax  be  lower  and  upper  bounds  on  the  velocity  c\  on  the 
interval  [0,D]:  for  all  2  G  [0,D], 


0  Cmin  —  £1(2)  5;  Cmax- 

We  will  require  that  8c  be  “sufficiently”  small  in  H 2:  We  choose  r  >  0  small  enough 
that  if  ||  c2-  Ci  ||H2[0)D]  <  r,  then 

(^)  II  ^2  —  <t>\  IIl°°[o,d]  <  2’ 

(11)  II  C2  ll//2[0,£)]  —  2  ||  Ci  ||J?2[0j£)] 

(12)  - - <  — -  <  2 -  for  all  0  <  2:  <  D,  and 

"  C"max  C2\Z)  C min 


(13)  |c2(2t)  —  ci(2)|  <  c2(z)  +  ci (2)  for  all  0  <  z  <  D. 

These  assumptions  assure  us  that  the  characteristics  t  =  (f>i(z)  corresponding  to  the 
two  velocities  are  not  too  far  apart  and  that  the  surface  measurements  depend  on 
comparable  depth  intervals. 


2.2.  Energy  identities. 
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2.2.1.  The  standard  energy  identity.  Given  a  velocity  c(z )  and  smooth  func¬ 
tions  u  and  v,  we  define  the  “standard”  energy  form  y  to  be 


V  = 


\dtu  dtv  +  dzu  dzv 
L  c2 


dz  +  [dtu  dzv  +  dzu  dtv]  dt. 


Then 


drj  =  —  ( dtu  □  v  +  u )  dz  A  dt. 

Now  suppose  r,  0  <  z0  <  z\,  and  T  =  r  +  2<f>(z1)  are  given,  and  define  T  to  be 

T  =  {{z,t)\ z0  <  z  <  zi,  r  +  (f>(z )  <t<T-  <f>(z)}  . 

We  define  the  operators  D+  and  D~  by  D±  =  ±(l/c)  dt+dz\  these  are  the  tangential 
derivatives  along  the  characteristics.  If  we  apply  Stokes’  theorem  to  rj  over  F,  then  we 
obtain 


(14) 


l 


ZO 

T-<f>(z0) 


J*1  (D+u  D+v)  (z,  t  +  <f>(z))  -  J*1  (. D~u  D~v)  (z,  T  -  <j>(z)) 


dt  [dtu  dzv  +  dzu  dtv ]  (z0,  t)  —  [  dz  f 

J  Zn  Jt 


T-0(z) 


Jt+4>(z0) 

In  the  case  u  =  v  we  obtain  the  standard  identity: 


r+4>{z ) 


dt  ( dtu  □  v  +  dtv  □  u) 


(15) 


JK  {d*u(z,t  +  <t>(z)))  -  j*  (n  u(z,T-  4>(z))y  = 

fT-<t>(zo)  rT-4>(z) 

/  dt2dtu(z0,t)dzu(z0,t)  -  /  dz  dt2dtuUu. 

Jt+4>(z0)  Jzn  Jt+6(z) 


2.2.2.  The  sideways  energy  identity.  Next,  for  a  smooth  function  u,  we  define 
the  “sideways”  energy  Qu  on  T  to  be 


rT—<t>(z)  . 

Qu{z)  =  /  dt  \(dtu(z,t))2  +  c2(z)  (dzii(z,t)): 

Jr+<h(z)  l 


fT-4>(z 
!t+4>(z) 

and  the  “sideways”  energy  form  u  to  be 

to  =  2  [ dtudzu \  dz  +  [(9<u)2  +  c2  (dzu)2 j  dt. 

Then 


du 


2 cc  (dzu)2  —  2 c2(dzu  □  u) 


dz  A  dt. 


If  we  apply  Stokes’  theorem  to  u  over  T,  then  we  obtain  the  sideways  energy  identity: 
(16)  Qu{zq)  =  c(z)  (d+u(z,t  +  <t>(z))Y  +  c(z)  (d~u(z,T  -  <t>(z))Y 


+ 


r d *  s 

J  zq  Jt 


T-4>(z) 

+<t>{z) 


dt 


2 cc'  ( dzuf  —  2 c2(dzu  □  u) 
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2.3.  The  progressing  wave  expansion.  In  this  section  we  will  discuss  the  pro¬ 
gressing  wave  expansion  for  solutions  of  the  wave  equation.  Associated  with  the  two 
velocities  c1  and  c2  are  the  two  wave  operators 


□  , 


dl  -  dl 


We  will  first  examine  the  solution  ut  of 


Ui  =  0,  z  >  0 
dzUi(0,t)  =  -6(t) 

Ui  =  0,  t  <  0. 

Recall  fa,  the  travel-time  to  depth  z  for  the  velocity  c,(z): 

«U  =  l‘d 

Jo  a(0 

For  smooth  c,-,  the  progressing  wave  expansion  [2]  shows  that  iq  will  have  a  jump 

along  the  characteristic  t  =  cj>i(z)  of  magnitude  a,(z)  =  cf(0)cf(z).  The  progressing 
wave  expansion  expresses  this  in  what  amounts  to  a  Taylor’s  series  expansion  valid  on 
one  side  of  the  characteristic  t  =  q f>i(z): 


N 


U 


;(M)  =  Y^Vk(z)Pk{t  -  (f>i(z))  +  RN+1(z,t). 

k= 0 


The  Pk  are  given  by 


Pk(x) 


0  if  x  <  0 
xk  if  x  >  0. 


The  coefficient  po(z)  is  the  jump  along  the  characteristic  t  = 

p0(z)  =  cf  (0)cf  (z), 

while  the  higher-order  coefficients  are  found  by  solving  the  second  transport  equation, 

(17)  2 p'ktfi  +  Pk4>"  ~  p'l- 1  =  0. 

We  can  easily  solve  (17).  From  the  condition  dzUi(0,t)  =  - 6(t )  we  find  that 

Pk(  0)  =  Ci(0K_r(0), 

and  so 

as)  Pk(z)  =  chupwpuw  +  ijf'chL, 

=  r/O’A  +  >PO)p’t  ,10)  !  -  (4  ) 
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The  remainder  i2/v+i(z,£)  is  the  solution  of 

Rn+i  =  Pn(z )  pN{t  ~  M2)),  z  >  0 

9zRn+i{0,  t)  =  p'N(0)PN(t) 

Rn+ i  =  0,  t  <  0. 

For  convenience,  we  will  denote  by  a,-  and  bi  the  first  two  coefficients  in  the  pro¬ 

gressing  wave  expansion  of  up. 

(19)  ai(z)  =  cf  (0)cf  (z) 

(20)  bi(z)  =  ^cf(z)  [cl  (z)a't(z)  +  cf  (0)a'(0))  -  jf  (cf)  o-  . 

One  interpretation  of  the  progressing  wave  expansion  which  will  prove  useful  is  that 
U{  is  the  solution  of  the  characteristic-initial- value  problem 

□i  Ui  =  0,  t  >  <f>i(z) 

(21)  dzUi(0,t)  =  0,  t>  0 

Ui(z,<j>i(z))  =  ai(z). 

Higher  derivatives  of  U{  solve  similar  problems;  we  will  also  use  the  fact  that  vt  =  dtUi 
is  the  solution  of  the  characteristic-initial- value  problem 


Vi 

=  o, 

t  >  <f>i(z) 

(22) 

dzVi(  0,  t) 

=  o, 

t  >  0 

Vi(z,4>i(z)) 

= 

while  Vi  =  dftU{  solves 

□iK- 

=  o, 

t  >  <f>i(z ) 

(23) 

dzVt(0  ,t) 

=  0, 

t  >  0 

Vi(z,<t>i(z))  =  hi(z), 


where  the  characteristic  value  h{(z)  is  the  third  coefficient  in  the  expansion  of  up. 

(24)  hi(z)  =  ^c](z)  (cf(z)6j(«)-cf(0)6j(0))  (c,71)  K  ■ 

A  similar  expansion  exists  for  6u{,  the  solution  of  the  linearized  problem  (7)  about 


8ui  =  -8<t>i{z)a,i(z)8(t  -  <f>i{z))  +  {8ai(z)  -  bi(z)8<f>i(z ))  H(t  -  <f>i(z))  +  R(z,  t), 
where  H  is  the  Heaviside  function,  and 


8<f>i(z)  =  -  f 
Jo 


is  the  first-order  change  in  the  travel-time.  The  coefficient  8ai  is 


8cii(z) 


1  6c(0) 


rW  +  |cf(0) 


8c(z) 

4( z) 


and  is  the  first-order  change  in  a,-.  For  the  remainder  of  the  paper  we  will  denote  8iii 
simply  by  8u. 
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2.4.  The  idea  of  the  proof.  Comparing  terms  from  the  progressing  wave  expan¬ 
sions  for  u2,  «i,  and  8u ,  we  see  that 

(25)  II  —  al  Ili2[0,£)]  —  C  II  C2  _  C1  lll,2[0,L»]  5 

(26)  II  &2  ~  IIl2[0,£)]  —  ^  II  C2  —  C1  Ili/2[0)D]  ' 

In  these  estimates  we  are  bounding  the  pointwise  values  of  ct-  and  its  derivatives  in 
terms  of  higher  Sobolev  norms.  We  also  have 

||  a2  —  ax  —  8a  ||L2[0vD]  <  C  ||  c2  —  cx  ||/^i [0>Jc»J  • 

Applying  the  standard  interpolation  result  for  Sobolev  norms  we  obtain 

(2^)  ||  «2  —  ax  —  <$a  ||L2[0)jD]  <  ^  II  c2  —  ci  IIl2[o,d]  II  c2  —  ci  II//2[0ijd]  • 

The  constant  C  depends  on  ||  cx  ||jyi[0D], 

mx  =  min  Ci(z), 

0<2<Z>  v  n 

D ,  and  r,  while  the  constant  I<  depends  on  ||  cx  mx,  £>,  and  r. 

On  the  one  hand,  these  estimates  show  that  the  velocity  perturbation  8c  causes  only 
a  small  difference,  commensurate  with  the  size  of  8c ,  in  the  values  of  the  ut-  along  the 
characteristics  t  =  (j>i{z).  Unfortunately,  the  characteristics  associated  with  the  veloci¬ 
ties  Ci  have  changed  as  well,  which,  as  described  in  the  introduction,  results  in  a  severe 
nonlinearity.  In  particular,  u2  —  ux  and  u2  —  ux  —  8u  do  not  solve  characteristic-initial- 
value  problems  but  instead  only  initial-boundary- value  problems  with  distribution  data. 
This  makes  a  direct  estimate  of  u2  —  ux  and  u2  —  u\  —  8u  inaccessible. 

We  will  instead  prove  the  estimates  in  Theorems  1.1  and  1.2  indirectly.  If  we  can 
show  that 

(26)  (f ’  F(c2)  -  F(ci)) L2[0<T]  <  A  ||  c2  -  ci  ||tfi[0iD]  ||  /  ||L2[0jT] 

for  all  /  6  Cp°(0,  A),  then  we  can  conclude  that 

||  F(c2)  —  F(ci)  ||£,2[0)T]  <  I<  ||  C2  —  CX  1 1  jy  1  [0,  JP]  ) 

which  ;r.  the  Lipschitz  estimate  of  Theorem  1.1. 

Similarly,  an  estimate  of  the  form 

(29)  |(/,  F(c2)  -  F{ci)  -  DF[Cl}(8c))L2[0T]\ 

—  II  IIl2[o,d]  II  IIh2[o,d]  +  B  ||  8c  ||^2]  ||  /  ||L2 

for  all  /  G  C£°(0,  T)  yields  the  estimate  in  Theorem  1.2. 

The  idea  of  the  proof,  then,  is  to  derive  bounds  of  the  form  (28)  and  (29).  To  this 
end,  suppose  that  w  solves 

□  i  w  =  0,  z  >  0 

dzw(0,  t)  =  f(t),  feC?(0,T) 
w  =  0,  t>T  —  <j>i(z). 
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If  we  apply  the  Stokes’  theorem  identity  (15)  with  u  =  Ui  and  v  =  w  over  the  region 
{(M)|0  <  ^  <  D,<f>i(z)  <  t  <  2 we  obtain  the  following  formulae  for 
(fi  F(ci)) L2[o,t]- 

(30)  [T  dt  f(t )  dtUl(0,  t)  =  [Ddz  D+w  (z,  M*))  4-  («1  (*,  Mz))) 

Jo  Jo  az 


(31)  [T  dt  /(f)  0tu2( 0,  t)  =  [Ddz  D+w  (z,  <f>2(z))  4  (u2  (z,  M*))) 

Jo  Jo  az 

rD  rT-fa  (z) 

-  dz  dtdtu2(z,  t)U2  w(z.t). 

Jo  J<fo(z) 

From  these  identities  we  see  that  we  will  need  to  estimate  the  difference  of  the  traces 
of  solutions  of  the  wave  equation  restricted  to  the  two  nearby  characteristics  t  =  (j>i(z), 
i  =  1,2.  This  is  the  subject  of  the  next  sections. 


2.5.  Trace  estimates.  In  this  section  we  present  estimates  for  the  traces  of  so¬ 
lutions  of  the  wave  equation  along  characteristics.  The  results  of  this  section  are  prop¬ 
agation  of  regularity  results  for  the  wave  equation,  with  explicit  dependence  on  the 
smoothness  of  the  velocity  coefficient.  These  trace  estimates  show  that  the  solutions  of 
the  wave  equation  are  a  great  deal  smoother  than  the  usual  trace  theorem  would  lead 
us  to  believe. 

These  trace  estimates  hinge  on  two  facts.  The  first  is  that  the  wave  operator  in 
one  spatial  dimension  2  is  hyperbolic  in  both  t  and  z,  a  consequence  of  which  is  the 
sideways  energy  identity  of  Section  2.2. 

The  second  point  is  that  we  can  almost  factor  the  wave  operator  into  derivatives 
along  its  characteristics.  For  i  —  1,2,  let 

d+  =  -dt  +  dz 

= - dt  +  dz 

C{ 

These  are  the  tangential  derivatives  along  the  curves  t  =  r  ±  <f>i(z).  Then 


(32) 

(33) 


D+D-  =  +  %dt 

cf 

DTDf  =  -a,  -  4a, . 

c- 


One  consequence  of  these  identities  is  the  commutator  relation 


(34)  DfD~  =  D-Df  + 

CHZ) 

Equations  (32)  and  (33)  indicate  that  we  can  factor  the  wave  operator,  up  to  a 
lower-order  term,  into  tangential  derivatives  along  the  forward  and  backward  charac¬ 
teristics.  This  fact  will  enable  us  to  derive  sharp  estimates  in  these  directions  for  traces 
of  solutions  of  the  wave  equation. 
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We  will  use  the  following  conventions  for  the  remainder  of  this  paper.  Recall  the 
depth  D  and  r  from  Section  2.1.  Given  a  velocity  c,-,  i  =  1,2,  let 

mi  =  min  cAz) 

o <z<D  v  ' 

Mi  =  max  Ci(z). 

o <z<D  v  ' 

Ci ,  Kt ,  and  Bt  will  represent  various  constants  whose  value  in  the  proofs  that  follow 
will  change,  even  from  line  to  line,  but  will  be  consistent  insofar  as 


Ci 

=  Ci  \ 

I<i 

A 

II 

( II  ci  II//2[0)dj  >  mi->  Mi,  D,  r\  , 

and 

Bi 

=  Bi 

(ll  ci  ll/f3[o,D]  i  mi,Mi,D,  rj  . 

The  dependence  on  the  maximum  Mi  is  mentioned  only  for  clarity  since  we  can  replace 
Mi  in  this  one-dimensional  problem  with  ||  c;  ||Hi[0  Dy 

We  will  also  often  suppress  the  domains  over  which  norms  are  computed.  Unless 
otherwise  noted,  the  norm  of  functions  of  2  alone  are  taken  over  the  interval  [0,-D], 
while  those  of  functions  of  t  only  are  taken  over  the  interval  [0,  T]. 

Finally,  recall  b  from  Section  2.1.  Given  r  >  —bj 2,  set  T,  =  r  +  2 <j>i(D)  and  define 
the  domain  T  to  be 

r  =  rr  =  {(m)I  o  <z<d,t  +  4>Az)  <t  <Ti  -  (f>Az)}. 

Associated  with  T  is  the  sideways  energy  Quy- 

QuAz)  =  [  [  ]  dt  \{dtu(z,  t))2  +  c2(z)  ( dzu(z ,  t ))2] 

J r+4>i(z )  L  J 

We  begin  with  a  general  proposition  relating  the  sideways  energy  QUiT  of  a  function 
u,  the  tangential  derivatives  Dfu  and  D~u,  and 

PROPOSITION  2.1.  Suppose  that  u  is  smooth  on  T.  Then  there  exists  a  constant 
C{  such  that 


(35) 

Dfu 

(z.t+aw)|!1I(0D]< 

Ci 

D~  u  (z,  T{  -  <f>i(z))  l2[q  d]  +  ||  dzu{ 0,  t) 

L2[o,Td  +  1  U  llz,2 (r)  ) 

(36) 

D~  u 

(z,Ti  —  4>i(z))  < 

V  ’  rn  11  L2[0,D]  - 

Ci 

Dt u(z,t  +  <t>i(z))  l2[o  d]  +  ||  dzu( 0,  t)  | 

lz/2[o,Ti]  +  II  u  llz/2(r)  5  and 

(37) 

Qu,t 

< 

Ci 

Df  u(z,T  +  (f>i(z))  L2[0  D]  +  ||  dzu( 0,  t )  1 

Il2[o,t,]  +  II  aiu  ilz,2(r)  • 
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Proof.  Let  z  G  [0,  D}.  If  we  integrate  the  sideways  energy  form  of  Section  2.2.2 
around  the  boundary  of  {((,t)\z  <  (  <  D,t  +  &(()  <  t  <  Ti  —  </>,(£)},  we  obtain 

QuA*)  =  facAO 0)f  +  facA0 (u-«(c,Tj  -  MO))2 

fD  fT,-MC)  r  - 

+  d(  dt  {2cfafa(()(dcu((fa)2-2c2(()(dcu((,t)niU((,t))  , 

Jz  Jr+<j>i(C)  L 

and  thus 


QuAz)  <  II  cWl°°  Dfu(z,r  +  <£,-(.z)) 


Z,2 


+ 


+ 


fD 

L  * 


1  +  2 


<(0 


Ci 


*(C) 


Qu,  r(C)  + 


D  u  (z,  Ti  —  <f>i(z)) 


□  t  U  T2  . 
L°°  *  L 


Z,2 


Applying  Gronwall’s  inequality,  we  arrive  at  the  sideways  estimate 


(38)  QuAz)  <  A  [  Dfu  (z,  r  +  &(*))  |ia  +  1  (+  Ti  -  M*)) 

where 

1 4(0 


I/2 


+  II  □;  u 


i  u  II Z,2 


A  =  He 


111,00 


+ 


II  L~ 


)  exp  ^ 


1  +  2 


c*(C) 


L1  ■ 


Meanwhile,  the  standard  identity  (15)  is 


(39) 


whence 


J^dz  (pfu{z,T  +  <t)i(z)))2  =  J^dz  [Di  u(z,Ti-  fi(z)))2 

fT  rD  tTi-4>i(z) 

/  dt  <9zu(0,  f)<9tu(0,  t)  +  dz  dt2dtu(z,t)Oiu(z,t). 

JO  JO  J T -\~4>i  [z) 


jQ  dz  (Dtu(z^{z)))2  < 

|  Dfu(z,Ti  -  <f>i(z))  L2  +  JQ  dt (dzu(0,t))2  +  a2  (dtu(0,t)f 


+ 


rD  rT-4>(z) 

dz 

Jo  J  <j>{z) 


dt 


r  l 


L  a1 


(□,  u(z,  t)f  +  a2  (d tu(z,  t))d 


where  a  is  a  real  number  whose  value  will  be  chosen  shortly.  From  this  inequality  we 
immediately  obtain 


2 

z,2 


jo  dz  (pf  u(z,T  +  (j>i(z)fj  2  <  |  Dtu  (z,  Tt  -  <j>i(z)) 

f  (ll  3M0,  t)  llh  +  II  □,  U II h)  +  Q2  (q,A o)  +  fa  <2„,rw)  • 
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Now  we  use  the  bound  on  Qu>T(z)  from  (38): 

(40)  JQdz  (Dfu(z,T  +  4>i(z))')2  < 

|  DJ u  (z,  Ti  -  <f>i(z))  fL2  +  (||  dzu( 0,  t )  ||22  +  ||  ||^) 

+  a2(l  +  D) A  Df  u(z,r  +  <j>i{z))  ||L2  +  ||  D~ u  (z,  Ti  -  <f>i(z))  |  ^  +  ||  u  \\2L2  . 

Set  o?  —  |  ((1  +  _D)A)_1.  Then  (40)  yields 

\j0dz  (Dtu(z,T  +  Mz)))2  <  \  |  D~ u  (z,  Ti  -  Mz))  fL2 
+  ^  (ll  d*< 0,  t)  W\2  +  ||  «  ||22)  +  i  II  u  \\\2 , 

from  which  (35)  follows. 

In  exactly  the  same  way  we  can  start  with  (39)  and  arrive  at  (36): 

Jq  dz  (p~ u  (z,  Ti  -  Mz))Y  < 

Ci  D?u(z,t  +  <f>i(z))  fL2  +  ||d*u(0,i)  11^2  +  ||  □,-u||i2  . 


Substituting  this  last  estimate  into  (38)  we  obtain 

Qu,r(z)  <  Ci  |  Df  u  (z,  t  +  </>,-(z))  ||^2  +  ||  dzii(0,  t)  ||22  +  ||  u  ||^2  , 
which  is  (37).  □ 

Recall  U{,  the  solution  of  (6)  with  c  =  C{.  We  will  use  the  preceding  proposition 
with  t  —  0  to  bound  the  energy  of  u,-,  dtUi ,  and  d2tUi  on  vertical  strips. 

PROPOSITION  2.2.  There  exist  constants  Ci,  Ki,  and  Bi  such  that  for  all  z  G  [0,  D), 


(41) 

(42) 

(43) 


f2<t>i(D)-4>i(z) 


r2<j>i(D)-4>i(z) 
r2  4>i(D)-4>i(z) 


dt  (dtui(z,t))‘2 


dt  (d2tUi(z,t)y 
dt  ( dzmUi(z,t))‘ 


<  Ci 

<  I<i 

<  Bi 


Proof.  In  Section  2.3  we  saw  that  Ui,  dtUi,  and  d2tU{  solve  the  characteristic-initial- 
value  problems  (21)  -  (23),  where  the  values  along  the  characteristic  t  —  f>i(z)  are 
explicitly  known.  This  together  with  (37)  yield  (41)  -  (43).  □ 

The  estimate  (41)  shows  that  the  map  F  :  ch  <9*u(0,  t )  extends  as  a  map  H1  — >  L2. 
Thus  the  boundary  value  <9tu(0,f)  in  problem  (6)  with  distribution  datum  is  defined  for 
velocities  in  H\ 
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The  estimates  of  the  Proposition  2.1  can  be  both  simplified  and  strengthened  for  a 
class  of  solutions  significant  to  our  proof. 

PROPOSITION  2.3.  Suppose  that  u  is  smooth  and,  that  u  vanishes  for  t  >  T  — 

Then  there  exists  a  constant  Ci  such  that 


(44) 


U  (z,  T  T  <f>i(z ))  5:  Ci  [j|  dzu(0,t)  ||L2[0)Ti]  +  II  171 1 


lL2(r) 


(45) 


L°°[0,£)] 


—  Ci  ||  dzu(0,t)  ||i2[0)7’i]  +  ||  Cfi  u  |  L2(r) 


(46) 


D •  u 

<  Ci 

% 

L°°(r)  —  2 

&U(M)  \\L2[0<Ti]  +  sup  j|  u(z,a  +  <j>i(z))  ||i2[0)jD] 

(J  >T 


Proof.  In  Section  2.1  we  arranged  matters  so  that 


Ti  -  Mz)  >T-  Mz) 


for  all  0  <  z  <  D,  so  D{  u  (z,  Tt-  —  <f>i(z ))  =  0.  The  estimates  (44)  and  (45)  then  follow 
from  (35)  and  (36). 

In  order  to  bound  the  T°°-norm  of  D~ u,  we  begin  with  (32).  For  (z,t)  E  T,  let 
t  =  a  +  <f>i(z).  From  (32)  we  have 


DfD-u  =  S dlU  (D}u  -  D-u)  -  O,  u, 

which  we  write  out  as  an  ordinary  differential  equation  along  a  characteristic: 


d_ 

dz 


_  I  1  cf. 

Pi  u{z,cr  +  Mz))\  +  2  f  u(z’a  +  Mz))  = 

1  d 

2  ~A+  uiz > a  +  Mz))  ~  ni  u(zi a  +  Mz))- 


l 

An  integrating  factor  for  this  differential  equation  in  D~u  is  cf : 


d_ 

dz 


’  i 

ci  (z)D~ u(z,cr  +  4>i(z )) 


(cf)'(z)Dfu(z,a  +  (f>i(z)) 


cf(z)aiU(z,a  +  </>i(z)). 


Since  the  depth  D  is  deep  enough  that  Di  u(D,a  +  <f>(D ))  =  0,  upon  integration  we 
obtain 

J(z)D~u(z, t)  =  JdC  (ci)'( QDfw  ((,  c  +  ^(C))  -  ci(0°i » «,  <r  +  MO) 

and  thus 


D :  u 


L°°(r) 


<  Ci  sup 


<7>T  L 


D?u(z,a  +  (j>i(z)) 


\l2  [o,r>] 


+ 1|  aiU(z,(T  +  <f>i(z)) 


\L2[0,D] 
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Then  (46)  follows  from  (44),  since  the  latter  assures  us  that 

sup  |  Dfu(z,a  +  4>i{z ))  \\l2[od]  <  Ci  [||  d zu{0,t)  ||i2[0Tt]  +  ||  □;  u  ||L2(r)]  . 

□ 

We  will  now  begin  to  apply  the  preceding  estimates.  Recall  w,  which  is  the  solution 


of 

□i  w  =  0 

(47) 

dzw(V,t)  =  m,  f  €  Co°(0,  T) 

w  =  0,  t>T 

-  Mz)- 

Proposition  2.4.  There  exist  constants  C\ 

and  K\  such  that 

(48) 

Df  w(z,t  +  4>i(z)) 

1  V  >  \  ))  l2[o,d] 

< 

Ci  II  /  II 

Il2[o,t] 

(49) 

Qw,t 

L°°[ 0,D] 

< 

Ci  II  /  II 

Il2[0,T] 

(50) 

DTw 

1  L°°(r) 

< 

Cl  II  /  II 

II2  [o,t] 

(51) 

\\dtw(z,T  +  Mz))\\L2[0,D] 

< 

Cl  II  /  II 

L2[0,T] 

(52) 

D7  D7w 

1  1  L~(r) 

< 

tflll/l 

Il2[o,t]  1 
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Then  from  (51), 


D^D-w  <7^11/11 


L2  5 


which  is  (52).  □ 

We  can  also  estimate  traces  along  curves  which  are  not  characteristics;  in  particular, 
we  can  obtain  trace  estimates  for  solutions  associated  with  C\  along  the  characteristics 
associated  with  c2,  and  vice  versa. 

PROPOSITION  2.5.  There  exists  a  constant  C\  such  that 


(53) 

(54) 


D%W  (z,  f2(z))  |L2j0)Dj  <  Cl  ||  /  ||L2[0iX] 
||  dtw  (z,  <f>2{z))  ||L2[0iD]  <  Ci  ||  /  ||L2[0)T] 


Proof.  We  will  prove  (53)  first.  Apply  Stokes’  theorem  to  the  standard  energy  form 
of  Section  2.2.1  over  the  region  {(z,t)\ 0  <  z  <  D,  <j>2(z)  <t<  <f>i{D)  +  <f>2(D)  -  <j>i(z)}. 
Since  <j>i(D)  +  <j>2(D)  -  <f>i{z)  >  T  -  <f>i(z)  for  0  <  z  <  D, 

Df  w(z ,  <f>i(D)  +  f2(D)  -  Mz))  =  0 

and  we  obtain 

4>i(D)+MD) 

dt  dzw(0,t)dtw(Q,t) 

=  1%  (d}w(z,m*)))2  +  (ww)J  -  f)  (d,w(z,U‘))f- 

Since  Df  —  Df  =  (<j>2  +  ff)  dt,  the  preceding  yields 
<h(D)+MD) 

dt  dzw(0,t)dtw(0,t) 

=  l  dz  (Dt™(z,<f>2(z))y  +  ^ ^  (D+w{z,<j>2(z))  ~  Dfw(z,<j>2(z)))2 

-  JQ  dz  (Dtw  (z,  Mz)))  -  +  ((-02  w  ( z ’  Mz)))  +  (jDfw(z,f2(z))y^ 

By  hypothesis  (13), 

#  -  ^2 

4>l  +  <t>2 

SO 

2  j0dz  {Dtw(z,Mz)))  <  jQdz  [Dfw{z,(j>2(z)))  +  Jo  dtdzw(0,t)dtw(0,t). 

The  estimate  (53)  then  follows  from  the  estimates  in  Proposition  2.4.  To  obtain  (54), 
note  that  D2  —  Dl  —  (<f>'2  +  <t>'i)dt  and  use  (43)  and  (53).  □ 
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Using  Propositions  2.4  and  2.5,  we  can  estimate  the  traces  of  solutions  in  which  w 
appears  as  a  datum. 

Proposition  2.6.  Suppose  that  v  solves 

□2  u  =  a(z)dtw 

(55)  dzv(0  ,<)  =  0 

»  =  0,  t  >  T  —  (f>\{z), 


where  a(z )  is  a  smooth  function. 

Then  there  exist  constants  C2  and  K2  such  that 


(56) 

D}v(z,  t  +  AM)Um 

< 

C2\\a 

IIl2[o,d] 

II  f  IIl2[o,t] 

(57) 

Qlr 

L°°[0,D] 

< 

C2  ||  a 

II  L*[0,D] 

II  f  Hl2[o,t] 

(58) 

Dfv 

2  l°°(  r) 

< 

K2\\a 

IIl2[0,D] 

II  f  Hl2[o,x] 

(59) 

II  (Z1  4*1  (z))  IIl2[0,£>] 

< 

C2  ||  a 

1  L2[0,D] 

II  f  Hl2[0,X] 

Proof.  The  bounds  (56)  —  (58)  follow  from  the  estimates  in  Proposition  2.3,  2.4, 
and  2.5.  The  key  is  that  the  bounds  in  Propositions  2.4  and  2.5  for  w  enable  us  to 
bound  the  A2-norm  of  dtw  and  thereby  □;  v  both  on  T  and  along  characteristics.  The 
proof  of  the  noncharacteristic  trace  estimate  (59)  follows  the  same  line  of  argument  as 
the  proof  of  Proposition  2.5.  □ 

Finally,  there  is  one  other  function  whose  traces  we  will  need  to  examine.  This 
function  arises  when  integrating  the  function  w  along  characteristics.  For  (z,£),  0  < 
z  <  D,  define 

(6°)  W(z,t)  =  d(  w  ((,t  -  ^(z)  +  . 

This  is  the  integral  of  w  along  the  characteristic  through  (z,t).  As  we  would  expect, 
it  is  one  degree  smoother  than  w  in  one  characteristic  direction  and  as  smooth  as  w  in 
the  other. 

PROPOSITION  2.7.  There  exist  constants  C\  and  K\  such  that 


(61) 

D!DtW(z,T  +  M*))\\L,M 

< 

^ll 

I/I 

Il2[o,t] 

(62) 

DfDfW 

1  1  1  L°°(r) 

< 

I<i\ 

1/ 

IL2[o,x] 

(63) 

DtWfrr  +  M*))  L(0,o] 

< 

Cl  II  / 1 

It2[o,x] 

(64) 

< 

Ah  | 

1  /  lln2[o,x] 

Proof.  The  bounds  (61)  and  (63)  follow  from 

D+D+W(z,t)  =  Dfw(z,t ) 
D+dtW{z,t)  =  dtw(z,t) 
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and  the  estimates  in  Proposition  2.4. 

Z^-bounds  on  D ^  W  and  D \  W  follow  from  the  fact  that  if 

U(z,t)=  f  d(u  ((,  t  —  <j>i(z)  +  ^i(C))  5 

J  D 

then 

Di U{z , t)  =  jT  dC  4(C)u  (C,  t  -  Mz)  +  MO) + cx(C)A"“  (C, t  -  Mz)  +  *(0)  • 

Finally,  the  bound  (64)  follows  from  the  identity 

D;d,w(z,t)  =  d;  (d*  -  n;)j  w(z,t) 

=  |  K»(M)  -  nrz>rw(*,i))  + 1  («,(»,  o  +  nfW(2,o) 

□ 

2.6.  Estimates  for  the  difference  of  traces.  The  estimates  in  this  section  are 
central  to  the  proof.  Here  we  will  derive  estimates  for  the  difference  of  solutions  of 
the  wave  equation  restricted  to  two  nearby  characteristics.  We  will  do  this  by  changing 
variables  from  (z,  t )  to  the  coordinates  corresponding  to  the  characteristics  t  =  T±<j>i(z). 
We  then  apply  the  trace  estimates  from  the  previous  section,  which  control  derivatives 
in  these  directions. 

To  this  end,  let 

£  =  II  4*2  —  4>\  IIl°°[0,D]  • 

For  0  <  z  <  D, 

\Mz)  -  Mz)\  =  [  —  <  —  Vd  ||  c2  -  d  \\L2 , 

Jo  C2  Ci  CiC2  L<x>  L 

so  invoking  assumption  (12)  we  obtain 

(65)  £  =  \\  4*2  —  4>l  Ilnoc^D]  <  Cl  ||  C2  -  Cl  || ■ 

Now  we  will  show  that  the  difference  of  the  traces  of  solutions  of  the  wave  equation 
along  the  two  nearby  characteristics  t  =  <f>i(z),  i  =  1,2,  is  also  small. 

Proposition  2.8.  Suppose  that  u  is  smooth.  Fori  =  1,2,  let 

Ai=  sup  (|  Dfu  (z,t  +  I  +  D+u(z,T  +  <l>i{z))  ) 

n*=  S^P,  (\\DtDtu(z,T  +  Mz))\\r2+  Df D~ u(z,r  +  4>i{z))  + 

|  D~ Df  u{z,r  +  <f>i(z))  \\l2  +  |  D~ D~ u{z,r  +  <f>{(z))  . 
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Then 


(66) 

II  «(*,&(*)) 

and 

(67) 

II  u(z,(f>2(z))  -  u 

u{zi4> i{z))  IIl2[o,d]  —  CiAi  ||  c2  ci  ||L2[0jD] 

(z,  ^(z))  -  (fa(z)  -  fa(z))  dtu  (z,  <fa(z))  ||L2[0>D] 

—  II  c2  C1  |Il2[0,D]  * 


Proof.  We  will  use  the  following  integral  form  of  Holder’s  inequality.  For  1  <  p  <  oo, 


(68) 


[//  u(z ,  t)  dt\  dz  1  P  <  J  [/  HM)|P  dz 


dt. 


For  simplicity,  we  will  prove  the  case  where  i  =  1;  the  case  i  =  2  is  the  same, 
mutatis  mutandis.  Let  fx(z)  =  f2(z )  —  <t>\ (z).  Then 


II  u(z,  Mz))  ~  u(z>  Mz))  II l2  = 


Cdz  [L 


n(z)+4>  i  (z) 


<t>l(z) 


dt  dtu(z ,  t ) 


fD  (  rn(z)+<h(z)  ci(z)  /  ,  x 

/o  *  (1W  *  2  “(M)  " 


Now  make  the  change  of  variable  from  t  to  the  r  coordinate  defined  by  the  charac¬ 
teristics:  t  =  r  +  <j>i(z).  If  we  do  this  and  apply  (68),  we  obtain 


u{z,Mz))  ~  «(*,&(*))  ||£a 


< 


JQdz  ^ fQ  dr  (Dtu(z,  t  +  <t>i(z))  —  D1  u(z,  r  +  ^>i(z)))^ 

II  Cl  IL~  II^IUoo  _sup  (|  Dfu(z,T  +  fa(z))  \\L2  +  |  D  +  u(z,T  +  <f>!(z))  ||L2)  . 


Then  from  (65), 


II  «(*,  Mz))  ~  u(z>  Mz ))  l|L2  <  C'iAi  ||  c2  -  d  ||L2 . 

The  proof  of  (67)  follows  exactly  the  same  path,  starting  with  the  identity 

=  j(A+n,l  -  A7 Dt  -  DfDT  +  D-D-), 

which  follows  from  dt  =  (ej If-  —  A7)-  We  then  proceed  as  above  in  order  to  bound 

“(*>&(*))  -u(z,^i(*))  -  (Mz)  ~Mz))dtu{z,<l>i{z))  =  /  dr  dsdlu(z,s). 
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Next  we  apply  the  previous  proposition  to  the  functions  considered  in  Section  2.5. 
PROPOSITION  2.9.  There  exist  constants  C\  and  K\  such  that 


(69) 

II  w  (z,  fa(z))  -  w  (z,  fa(z))  ||L2[0(£>] 

< 

C-il 

|  C2  -  Cl 

IIl2[o,d]  II  f  IIl2[o,t] 

(70) 

I  D1  w  (z,  ^(z))  -  Df  W  (z,  fa(z))  L2[0)D] 

< 

K\ 

||  C2  -  Cl 

IIl2[o,t>]  II  f  IIl2[o,t] 

(71) 

\\dtW(z,fa(z))-dtW(z,fa(z))\\L2[’M 

< 

Ci| 

I  c2  -  Cl 

IIl2[o,d]  11  f  IIl2[o,t] 

and 


II  W  (z,  <f>2(z))  -  W  (z,  fa (z))  -  (fa(z)  -  fa (z))  dtW  (z,  fa (z))  \\L2[0>D] 
<  K\  ||  c2  —  ci  ||L2[0iij]  ||  /  IIl2 [o,t]  • 

If  v  solves 


□2  v  =  a(z)dtw 


(73)  <92u(0  ,f)  =  0 

v  =  0,  t  >  T  —  fa (z), 

where  a(z)  is  a  smooth  function ,  then  there  exists  a  constant  K2  such  that 

(74)  ||  v  ( z ,  fa(z))  -  v  (z,  fa (z))  ||L2[0  D]  —  7f2  ||  a  ||L2[0iD]  ||  /  1 1 £,2 [o,x]  ■ 


Proof.  We  apply  Proposition  2.8  and  the  trace  estimates  of  Section  2.5. 

To  obtain  (69),  we  use  the  estimates  on  Dt  w  and  Dr  w  in  (48)  and  (50).  Estimate 
(70)  follows  from  (54)  and  (52),  since  Df  Dfw  =  (c[/cl)dtw  (see  (32)).  Similarly,  (71) 
follows  from  (61)  and  (62).  The  second-order  estimate  (72)  follows  from  Proposition 
2.7. 

Finally  (74)  follows  from  (56)  and  (58). 

□ 

2.7.  A  bound  on  area  integrals.  Besides  estimating  the  difference  of  traces 
along  characteristics,  we  will  also  need  to  estimate  differences  of  the  form 


rD  /•T-0i(z) 

dz  dta 

Jo  J<t>2(z) 


rD  rT—  4>\{z) 

(z)dtwdftu2—  dz  dt  a(z)  dtw  d?tui. 

Jo  J<t>\(z ) 


Here  we  will  see  the  first  application  of  the  estimates  for  the  difference  of  traces  from 
Section  2.6.  We  have  the  following  estimate. 

PROPOSITION  2.10.  There  exists  a  constant  B\  such  that  for  any  smooth  a  =  a(z), 


(75) 


rD  rT—<t>\  (z)  rD  rT-<h(z) 

dz  dt  a(z)  dtw  d:tu2  —  dz  dt  a(z)  dtw  dfux 

Jo  JMz)  Jo 


<  B 1 1|  c2  -  ci  ||H, 


lL2 


01  (z) 

f  I 


lL2 


R.  M.  Lewis  and  W.  W.  Symes  -  Solutions  of  the  wave  equation 


24 


Pr:of.  Let  v  be  the  solution  of 

□2  v  =  a(z)dtw 
dzv(0 ,t)  =  0 

v  =  0,  t  >  T  —  <j>i(z). 

Now  we  apply  the  Stokes’  theorem  identity  (15)  to  the  two  regions 
{(z,t)\0  <  ^  <  D,<f>i{z )  <  t  <  2 (f>i(D)  -  <f>i(z)} 
for  i  —  1,2.  In  the  case  i  =  2  we  have 

fD  rT—4>\(z)  i-D 

dz  dt  a(z)dtw  d^tu2  =  —  /  dz  D%v  (z,  <j>2(z))  Dtu2  (z,  Mz)) . 

Jo  Jo 

On  the  other  hand,  since 


□2  w  =  (m2  —  Qi  )tw  = 


5  -  3  * 


the  Stokes’  theorem  identity  (15)  in  the  case  i  =  1  yields 


D,  i(2)  „  fu  r-i-<Pi(z) 

'z  /  dt  a(z)dtwdftui  =  dz  dt  ((n2  -  Dj )  +  Dj )  v  dftui 

j<t>l(z)  Jo  J<t>i(z) 

fD  fT-4>  i(z)  /  1  \\  n  rD 

=  Ldz  L  I  \  dt  \  72  ”  3  )  dttv  dttu  1  -  dzD+v  (z,  (f>x(z))  D+  <9tUi  (2,  ^a(ar)) . 

*/0i  (2)  \^2  Ci  y  ./o 


j  dz] 


rD  fT—<j>i(z) 


Recall  from  the  progressing  wave  expansion  that  d?tUi  (z,  <f>\{z))  =  bi(z)  (see  (22)). 
Then  we  may  integrate  by  parts  in  t  to  obtain 


fD  fT-4>i(z)  (  \  1 

dz  dt  —z - r 

Jo  J<t>i{z)  \C2  c{ 


r%  U- 


dlvdlu , 


dtv(z,(f>1(z))  61  (z,<f>x(z)) 


Jo  \4(z)  c^)/  ‘  v  ,riv  "  1V 

fD  fT-<fn(z)  (  \  1  \ 

-  dz  dt  [- 2  -  ~2  dtvdfttu 1. 

io  -/*(*)  Vc2  ci/ 


From  the  formula  (20)  for  &i  and  the  trace  estimate  (59)  for  dtv  (z,  4>i{z))  from  Propo¬ 
sition  2.6,  we  obtain  the  bound 


<  Ddz  (~ 


c\(z) 


-  dt^(z,Mz))bi(z,Mz))  <  I<\  ||  c2  -  Cl  ||L2  ||  a  \\L2  11/ 


Meanwhile,  from  the  sideways  estimates  (42)  and  (57)  we  obtain 


fD  fT-M*)i  (l  1 

dz  dt  [  —z - 7 

Jo  \C2  C\ 


cl  c 


^  dtv  dtttu!  <  Bi  II  c2  -  Cl  ||  jj2  II  a  \\L2  ||  /  ||L2 . 
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Thus 


D  fT-M z)  (  l 


f  dz  I 

JO  J<b 


dt 


<t>i  (2) 


9  ®ttv  ^ttu l 


<  Bi  c2  —  ci 


I L2  II  a  liL2 


L2 


Meanwhile,  since  n,-  (z,  <j>i(z))  = 

l-D  rD 

J ‘o  dz  v  (z,<f>2(z))  D^u-2  (z,  4>2 (z))  -  J ^  dzD^v(z,(f>1(z))  (z^^z)) 


=  /  dza^-Djvfo^z))-  /  dzaitzJDfvfo^z)) 

do  do 

-  J0dz  (a2  _  al)/  (2  W (2,  ^2(2))  ~  fQdz  ai(Z)  b  (Z>  ^2(2))  -  u  (z,  ^l(z))] 

Applying  (56)  to  the  first  integral  and  (74)  to  the  second,  we  obtain 

')tv(z^(f>2(z))  D^u2{z,4>2(z))  -  f  dz  D^v(z,</)1(z))  Dfui  (2,^(2)) 


<  Kx  II  c2  -  Ci 


\m  II  a  II l2 


L2  • 


3.  The  proof  of  Theorem  1.1:  Lipschitz  continuity.  We  will  now  prove  Theo¬ 
rem  1.1.  The  progressing  wave  expansion  tells  us  that  along  the  characteristics  t  =  <f>i(z), 

i  —  1,2, 


Ui(z,<f>i(z))  =  a,i(z). 

Thus  (30)  and  (31)  become 

/  dtf(t)F(c1)  =  j  dza[(z)Dfw(z,(j)1(z )) 
do  do 


and 


[l  rU  nD  rl  —<pi(z) 

I  dtf(t)F(c2)=  /  dz  a'2(z)D2W  (z,  4>2{z))  +  dz  dt  dtu2  Q2  to. 

v/o  ./ o  J  4*2(2) 

Once  again  we  rely  upon  the  fact  that 


□2  w  =  (D2  —□!)«;=  f  - 


c2  c2 
c2  ci 


We  use  this  and  integrate  by  parts  in  t.  Bearing  in  mind  that  w  =  0  for  t  >  T  -  4>\(z) 
and  d tu2  (z,<f>2(z))  =  b2(z),  we  obtain 

rD  rT—i)>i(z)  rD  rT-<t>i(z)  (  \  \  \ 

dz  dt  dtu2  d2  w  =  I  dz  I  dt  j  ^ - 5-  dlw  dtu2 

do  J<t>2(z)  Jo  VC2  ClJ 

-L 


D 

dz 


i(*) 


1  \  tD  rT—<t>i(z)  {  1  \  \ 

b2(z)dtw(z,<f)2(z))  -  dz  dt  [  — - -  j  dtw  d2ttu2 

z\\z) )  do  d02(z)  \c2  cl) 


R.  M.  Lewis  and  W.  W.  Symes  -  Solutions  of  the  wave  equation 


26 


and  so 

rT 


[  dtf(t)F(c2)=  [  dza2(z)D%w(z,<f>2(z)) 

Jo  Jo 

fD  f  1  1  \  [D  rT-M*) 

/  dz  ~YT\  -  TTT  b2{z)dtw(z,<f>2(z))  -  dz  dt 

Jo  \CZ2(z)  cf(z)J  Jo  Jfa(z) 


+ 


We  also  have 

rD 


A  c\) 


\ 


dtwdlu2. 


Jru  rjj 

dza'2(z)D$w(z,<f>2(z))  -  /  dz  a[(z)Df  w  (z,  <J>i(z))  = 

0  Jo 

[  dz  [a2(z)~a,1(z))]Dfw(z,4>2(z))-  j  dz  a”(z)[w  (z,<f>2(z))  -  w  (z,  <f>\(z))] . 
Jo  Jo 


Thus 

rT 


(76)  /  dtf(t)  [f(c2)  -  F (c)]  = 

Jo 

j  dz  [a2(z)  -  a[(z))]D^w(z,(j)2(z))  -  f  dz  a"(z)[w  (z,  (f>2(z))  -  w  (2,  </>i(z))} 

Jo  Jo 

fD  (  1  1  \  fD  fT—(f>i  (z)  /  l  l  \ 

+  Jo  dz\~c W)~W)}  ^z)diw{z’ -  Jo dz  JMx)  dt  U  -  4 )  diw9«u 2' 


We  can  bound  the  terms  in  (76).  From  Proposition  2.5  we  have  the  following  bound 
on  the  first  term: 

/  dz  [a'2(z)  -  a[(z))]  D^w  (z,<f>2(z)) 

Jo 

Using  (69),  we  bound  the  second  term  in  (76)  as 


<  Ci  ||  c2  —  Ci 


\m 


f 


\L2  • 


/  dza"(z)[w(z,<t>2(z))  -  w  (z,<J>i(z))] 

I  •J  0 


<  AT  ||  c2  -  Cl 


\l2 


f 


I L2 


The  third  term  on  the  right-hand  side  of  (76)  is  bounded  using  the  expression  for  b2  in 
(20)  and  the  trace  estimate  (54)  for  dtw  in  Proposition  2.4: 


Jo 


D 

dz 


1 


b2(z)dtw(z,<j>2(z)) 


<  I<2  ||  c2  —  ci 


f 


c2^  c2^  J  -1  ^2  ||  '-2  -  w  ||L2  ||  j  ||L2  • 

Finally,  we  estimate  the  area  integral  by  means  of  the  bounds  on  the  sideways 
energy  of  dtu2  and  dtw  given  in  (42)  and  (49): 


fdz  [ 


T-4>  i(z) 


dt 


$2(2) 


~2  -  4  )  dtw  d2ttu2 


<  K2  ||  c2  —  ci 


\L2 


f  I 


I L2  ■ 


Since  we  assume  that  ||  c2  —  Ci  ||ff2  <  r,  we  may  replace  the  constants  K2  with 
constants  AT  in  the  preceding  bounds.  We  conclude  that 

rT 


/  dtf(t)  [F(c2)  -  F(ci)\ 
Jo 


—  [-^1  II  C2  ~ '  C1  r,2  T  Cl  c2  —  Cl 


\m\ 


11/ 


\l2 


for  any  f  £  C£°[0,  T],  from  which  it  follows  that 

II  F(c2)  —  A(ci)  ||L2  <  [Ki  ||  c2  —  ci  ||L2  +  Ci  ||  c2  —  ci 
This  proves  Theorem  1.1. 


1  ml 
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4.  The  proof  of  Theorem  1.2:  Regularity. 

4.1.  A  formula  for  DF[c^\{8c).  We  begin  with  a  formula  for  DF[c\]{8c). 
Proposition  4.1. 

(77)  {/,  DF[c i](6c))i2[0T]  =  dz8a'(z)D+w(z,cj)i(z)) 

+  J  dz  a'[(z)8(j)(z )  +  bi(z)^j  dtw  (z,  <f>t (z)) 

rD  rT—4>\(z)  Of.  r(  z\ 

+  dz  dt  -  ~-dtw(z,t)d^iu1(z,t). 

Jo  J(j>  1(2)  c^yzj 

Proof.  We  will  need  the  following  integration  by  parts  identity.  Suppose  6  is  a 
smooth  function  such  that  6(0)  =  0.  Then  we  claim  that 

(78)  f  dz  (cl  b  \  (z)Df  w  (z,  <f>i(z))  =  f  dz  c^  (z)b'(z)-^—dtw  (z,<t>i{z)) 

Jo  \  /  Jo  cl[z) 

Identity  (78)  is  proved  as  follows.  From  (32)  we  have 

JQ  dz  c\  (z)b(z)  dtw{z,<f>1{z))-^D^w{<z,<j)l{z))  =0. 

Integrating  the  second  term  produces 

-  J0dz  cUz)b(z)-fDiw  (z,<h(z))  = 

L°dz  (^',(2)+c'(z)i,'(2)) 

Thus, 

JQ  dz  cf  (z)b(z)  C^~dtw{z,<f1{z))  -  -^-Dfw(z,<f>i(z))  = 

jQ  dz  (cl'Sj(z)b(z)Dfw(z,4>1(z))  +  dzcl(z)b\z)Dfw(z,<f>1(z ))  =  0, 

and  so 

Jo  dz  (cf'j  [z)b{z)Dfw{z,(j)1{z))  +  dz  c{{z)b'  [z)Df  w  [z,(j>i{z)) 

-  [  dzcl(z)b\z)—^—dtw(z,<j>i(z))  =  0, 

Jo  Ci yz) 

which  is  (78). 

We  now  turn  to  deriving  (77).  Applying  the  Stokes’  theorem  identity  (15)  and  a 
straightforward  limiting  argument  to  handle  the  singular  nature  of  the  integrand,  we 
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obtain 

(79)  J  dt  f(t)  dt8u(0,t)  =  J  dz  (6a  —  biS<j)),(z)Df  w  (z, 

+  (Ddz  fT~Mz) dt?^dtw(z,t)%Mz,t) 

Jo  J<h(z)  C{(z) 

We  then  apply  (78)  to  obtain  the  following: 

Jq  dz  (bi8<j))'  (z)D^w(z,<t>i(z))  =  jf  dzcl(z)  2b16<f?J  (z)-^j^dtw  (z,  (j>i(z)) . 
However,  we  see  from  the  formula  (20)  for  bi  that 


ci  2&i  )  (*)  =  \ci  2W<(^)> 


so 


Sc(z )' 


c\(z )  (cj  2M<£)  (z)  =  c\(z )  ^  2(z)a"(z)6(f>(z)-c12(z)-^J  b^z). 


Thus 

rD 


Jo  dz  (biStff)’  (z)Dfw  (z,<f>x(z))  =  JQdz  [o!{(z)H(z)  ~  dtw(z,<f>1(z)) . 


Proposition  4.1  then  follows  from  this  latter  equality  and  (79).  □ 

4.2.  Continuity  of  DF.  Proposition  4.1  shows  that  DF[c]  :  H1  — >  L2  for  c  G  H2, 
since  we  can  easily  use  the  trace  estimates  in  Proposition  2.4  to  derive  a  bound 


(/,  DF(c](6c)) L2[ qj’] 


<  ATI  8c  | 


m 


f 


L2 


for  all  /  €  Co°[0,  T],  where  K  —  K  (||  c  ||#2[0  d]  ’  cm»n,  D)  ■  However,  in  order  to  show 
that  DF[c]  depends  continuously  on  c,  we  must  require  more  smoothness  of  both  the 
velocity  and  the  admissible  perturbations  8c.  The  following  is  part  (ii)  of  Theorem  1.2. 
PROPOSITION  4.2.  There  exists  a  constant  B\  such  that  if  ||  c2  —  C\  ||#2[0£)]  <  r> 

then 


||  DF[c2 ]  -  DF[ci]  ||  <  Bx  j|  c2  -  Cl  \\H2[0>D] , 

where  the  operator  norm  is  taken  over  operators  from  H2[0,D]  to  L2[0,T]. 

Proof.  We  begin  with  the  identity  (79)  from  the  proof  of  Proposition  4.1:  For 

*  =  1,2, 


Jo  dt  f(t)dt8ui(0,t)  =  J  dz  (8ai  —  biSfi)' (z)Df  w  (z,<j>i(z)) 


+ 


rD  rT—<j) 

/  dz 

Jo  J<t>i(z ) 


T-M.)  .  2Sc(z) 


dt 


c?(z) 


■dtw(z,t)d„Ui(z,t). 
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Then 


(/,  (DF[c,}  -  UF[c])  (ic))^  = 

J  dz  ( 8a2  —  b2S(j)2),(z)D2W  (z,  <f>2(z))  —  J  dz  (8ax  —  b16<j)1)'(z)D^w  (z,  <j>i{z)) 

[D  fT-<h(z)  26c  _  _o  (D,  fT-Mz)  2 Sc.  -2 

+  dz  dt  -^-otw  cLu2  —  /  dz  /  at —z-otw  otu-y. 

JO  J<j>  2(2)  C2  JO  Cl 

.After  adding  and  subtracting  a  term  to  the  first  difference  we  can  integrate  to 
obtain 

J  dz  (8a2  —  b2S  fa)' (z) w  (z,  <j>i(z))  —  J  dz  ( Sai  —  b\8 <f>i)\z) D+ w  (z,  <j>i(z)) 

=  J  dz  [(<5a2  -  b28<j)2)\z)  -  (8ai  -  w  (z,  <f>i{z)) 

-  f  dz  (8^  -  bi6<f>i)"(z)  [w  (z,  <^(2))  -  (z,  <j>i{z))] . 

J  0 

Applying  the  trace  estimates  in  Propositions  2.4  and  2.9  we  can  bound  the  preceding: 

dz  ( Sa2  -  b2S<f>2)\z)D. w  (z,  <f>2(z))  -  dz  (8ar  -  b18<j>i)'(z)D f  w  (z,  <j>\(z)) 


<  K 1 II  c2  -  ci 


8c  | 


/! 


Il2 


Meanwhile,  from  Proposition  2.10  we  obtain  the  estimate 

rD  rT-fa(z)  2  8c  -  rD 


f  dz  [  iZ  dt  ~dtwd2ttU2-  f  dz  f 
Jo  J<t>2{z)  ci  Jo  j<t> 

■D 


< 


t-4>i(z)  28c  2 

dt  —j-otw  OttUi 


/*/ 

Jo  j<t> 

fdz  / 
Jo  Jd> 


T—<t>i(z)  28c 

dt  — — ' 

<t>2(z)  c2 


<h(z) 

D 


T-4>  iW 


:  „  ru  rz-viyz)  Jhr 

-d,w  d2uu2  —  dz  dt  -r-dtw  d2tu1 

Jo  J<h{z)  c% 

(2 8c  2 8c\  2 

cu  I  - 3-  I  Otw  OttUi 

<Br  ||  c2  —  Ci  \\L2  ||  8c  ||L2  ||/||L2. 

Thus  we  conclude  that 

(/*  (■QP’fcs]  ~  DF[c  1])  (<5c))L2[o:r]  <  Bx  ||  c2  -  ci  || H2  ||  8c  \\H2  ||  / 
for  all  /  E  C£°[0,  T],  and  therefore 

||  DF[c2]  -  i3F[ci]  ||  <  Bx  ||  c2  -  ci  ||„2  . 


z,2 
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4.3.  Estimating  the  linearization  error.  From  (76)  and  (77)  we  obtain 

(80)  </,  F(c2)  -  F(c)  -  DF[Cl}(6c))L2[0T]  = 
fD 

/  dz  (a2  —  ax  —  8a)' {z)D^w  (z^(j)2{z)) 

J  o 

fD  r  , 

+  Jo  dz6a'(z)  [. Djw(z,<j>2(z ))  -  («, 

-  f  dz  a"(z)  [w(z,(j)2(z))  -w(z,<f>1(z))  -  ( <f>2(z )  -  ^(z))  dtw  (z,  <j)\{z))\ 

J  0 

-  I  dz  a"(z)  [4>2{z)  —  4>i{z)  —  8cj)(z)]  dtw  (z,  <j>i(z)) 

J  0 

-  f  dz  {-y-T - y-T  +  2  3^^)  b2{z)dtw(z,<f>2(z)) 

Jo  \4(z)  cl(z)  c\(z)  ) 

+  f0  dz  (M*)  -  M*))  dtw  (z,  Mz)) 

2 Sc(z) 


L  dz  4(z)bl ^  ^tW  ^ Z ^  _  dtW  (Z’ 

&  Cw*  -  4 + a«“2 


rU  fl-<h  (z)  2dC  „  rD  Of,r 

+  dz  dt  ~dtw  dlu2  -  dz  dt  ~dtw  d2ttUl 

Jo  J<j>2(z)  c\  Jo  J<t>i(z)  c\  “ 

Estimates  for  most  of  these  terms  are  immediate.  From  (27)  and  (53), 


J  dz  (a2  —  ai  —  ba)\z)D2W  (z,  <f>2(z)) 


<  Cl  ||  C2  —  Cj  11^2  II  C2  —  Cl 


\m  ii  J  iil2 


Integration  by  parts  and  the  trace  difference  estimate  (69)  yields 


fQ  dz  ba\z)  [Dt w  ( z ,  Mz))  -  D+w  ( z ,  fa (z))] 
From  the  trace  estimate  (51), 

rD 


—  Fi  ||  c2  ci  W12  ||  c2  Ci  ||^i  ||  /  \\L2 


Jo  dza'l(z)  [(j)2{z )  -  <f>i(z)  -  84>(z)]dtw(z,<j>i(z)) 


<  I<1  ||  C2  -  Cl  \\\2  ||/||L2  . 


The  next  term  in  (81)  is  bounded  using  the  definition  of  b2  in  (20)  and  (54): 


r^J  . 1  . 


l Jo  'TTF  ■  W) +  dFT )  h{z)d'w' {z' Hz)) 

Using  (26)  and  (54),  we  have 
rD  2 8c{z) 


<  Ki  ||  c2  -  ci 


lL2  II  J  Iil2 


L dz  ( b ~ bi^ dtw  &  ^ z ^ 


<Al||C2  C\  \\L2  ||  C2  C\  ||^2  II  /  11^2 
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Finally,  from  (41)  and  (49),  which  bound  the  energy  of  dtu2  and  w  on  vertical  strips, 
we  have 


D  rT—4>i(z) 


l  dZf 

JO  J<t> 


dt 


<Mz) 


4  -  4  +  ^)  a,u,alm 


r2  r2  r3  I 

c2  cl  ci  / 


<  AT  II  c2  —  ci 


I L2 


L2 


Three  terms  remain.  The  first  of  these  is 

rD 


(81)  /  dz  a”(z )  [tu  (z,  fa(z))  -  w  ( z ,  Mz))  ~  (Mz)  ~  Mz ))  9tw  ( z ,  fa(z))\ 

J  0 

Recall  the  function  W(z,t)  from  (60)  and  note  that 

-^W  (z,fa(z))  =  w(z,t) 

-^-W  (z,fa(z))  =  (fa  -  fa)' (z)dtW  (z,fa(z)) . 

This  enables  us  to  integrate  (81)  by  parts: 

JQ  dz  ai(z)  lw  (zi  Mz))  ~  w  (z,  fa(z))  -  (fa(z)  ~  fa(z))  dtw  (z,  fa(z))\ 

=  J^dza"(z)  —  (W(z,fa(z)))  -  (fa  -  fa)' (z)dtW  (z,fa(z)) 

-  J^dza"(z)-^-(W  (z,fa(z)))  -  Jodza"(z)6(/)(z)^-(dtW  (z,fa(z))) 

=  -  /  dza'f(z)  [W  (z,fa(z))  -  W  (z,fa(z))  -  (fa(z)  -  fa(z))dtW  (z,fa(z))\ 

JO 

+  f  dz  a"(z)  (fa  6<f> )  (z)  dtW  (z,  fa (z)) 

J  0 

-  f  dz  a”(z)  (fa  -  fa-  8(f>)' (z)  dtW  (z,  fa  (z)) 

J  0 

+  [Ddz  a'l(z)  8<t>(z)  [dtw  (z,  fa(z))  -  dtW  (z,  fa(z))\ . 

J  0 


We  can  estimate  each  term  in  this  latter  expression.  By  the  second-order  trace  difference 

estimate  (72), 

/  dz  ci'"(z)  [W  (z,  fa (2))  -  W  (z,  fa(z))  -  (fa (2)  -  fa (z))  dtW  (z,  fa (z))] 

J  0 


<K  1 


(cir 


L2 


c2  Ci  ||L2 


L2  ■ 


Next,  the  definition  of  W  and  the  bound  on  dtw  (z,  4>i(z))  in  Proposition  2.4  enable  us 
to  bound  dtW  (z,  ^i(z)),  whence 


I  dz  a”(z)  (fa  -  fa-  8<f>)' (z)  dtW  (z,  fa  (2)) 

J  0 


<  AT  ||  c2  —  ci 


1 1? 


L2  ■ 
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Finally,  from  (71)  we  obtain 

/  dz  a"(z)8<j)'(z)  [dtW  ( z ,  <fo(z))  -  dtW  (z,  <t>i{z))\ 

J  0 


<  K\  II  c2  —  ci 


L*  II  j  III ?  ■ 


Thus, 


f  dz  a"(z)  \w  (Z,MZ))  ~W(Z,MZ))  ~  8(f>(z)dtw  (z,^(z))] 

J  0 

<  Bl  II  C2  -  C!  \\l  ||  /  ||„  . 

Next  we  will  estimate 

rD  2  6c(  z} 

/  dz  ~o7\^\{z)  [dfw(z,<t>2(z))  -  dtw(z,4>i(z))]  . 

J  0  (  Z  ) 


This  is  straightforward.  Since 


dt  =  -  (b+  -  Dr) 


CjC2 

Cl  +  c2 


{pi  -  d-,) 


we  have 


J0  dz  \.dtw  (z,  <j)2(z ))  -  d tw  (z,  <f>\(z))) 


c\{z) 

=  J0  dz  2'c3^jMz)^Y^  [Dtw(z,<j>2(z))  -  D~l w  (z,  </>i(z))] 

,  [Dj  2 8c(z)h  (  cx(z)cv{z)  ci(z)\  +  , 

+ Jo dz  ~^ybi{z)  w(7)+c2(z)  mz)) 

rD  2  8c(z)  ci(z) 

-  J0  dz  C?(z)  bl^~2  W  ~  D'  W  ^(7) 

fDj  2 8c(z)ur  s(  ci (z)c2(z)  ci(z)^  /  J  f 

Jo  dZ  c?(z)  bl ^  (ci(z)  +  c2{z)  2  J  ^  W  ' 

For  the  first  term  in  this  latter  identity  we  integrate  by  parts  and  apply  (69).  The  third 
term  we  bound  using  (70);  the  other  two  terms  are  bounded  by  the  trace  estimates  in 
Propositions  2.4  and  2.5.  We  arrive  at 


rD  2 8c(z) 

Jo  dz  ~4(z)  i bl^  ^dtW  ~  dtW  ^7))} 

<  I<  1  ||  C2  -  Cl  ||L2  ||  C2  -  Cl  ||H1  ||  /  ||La  . 

Finally,  applying  Proposition  2.10  to  the  difference  of  the  area  integrals  we  obtain 
the  bound 

D  rT-4>i{z)  28c 

dz  /  dt  —^Otw  ottui 

<  Bl  II  <2  -  Cl  II  /  ||i2 . 


fD.  [T-M*)  28c  r 

dz  dt  -^-dtw  dftu2  - 

Jo  Jfa(z)  C{  Jo 
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The  conclusion  is  that 


< 


{/,  Fid)  ~  C(c.)  -  flFkKfc))^ 

.Bl  II  C2  -  c,  Up  +  Ki  ||  Cj  -  d  Up  II  c2  -  c,  ||„,  I  II  / 


I L2 


for  all  /  G  Co°[0,  T].  This  estimate,  the  continuity  of  F,  and  the  continuity  of  DF  yield 
Theorem  1.2. 
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